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ABSTRACT 



< 
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' Context. Central starbursts in galaxies are an extreme example of ongoing galaxy evolution. The outer parts of galaxies 

contain a fossil record of galaxy formation and evolution processes in the more distant past. The characterization of 
resolved stellar populations allows one a detailed study of these topics. 

Aims. We observed the outer parts of NGC 1569 and NGC 4449, two of the closest and strongest dwarf starburst 
^ ' galaxies in the local universe, to characterize their stellar density and populations, and obtain new insights into the 

, structure, formation, and evolution of starburst galaxies and galaxy halos. 

, Methods. We obtained HST/WFPC2 images between 5 and 8 scale radii from the center, along the intermediate and 

' minor axes. We performed point-source photometry to determine color magnitude diagrams of / vs. V — I. We compared 

the results at different radii, including also our prior HST/ACS results for more centrally located fields. 
Results. We detect stars in the RGB and TP-AGB (carbon star) phases in all outer fields, but not younger stars such as 
those present at smaller radii. The RGB star density profile is well fit by either a de Vaucouleurs profile or a power-law 
^ ' profile, but has more stars at large radii than a single exponential. To within the uncertainties, there are no radial 

gradients in the RGB color or carbon-to- RGB-star ratio at large radii. 
0^ I Conclusions. The galaxies have faint outer stellar envelopes that are not tidally truncated within the range of radii 

00 , addressed by our study. The density profiles suggest that these are not outward extensions of the inner disks, but are 

■ instead distinct stellar halos. This agrees with other work on galaxies of similar morphology. The presence of such halos 

' is consistent with predictions of hierarchical galaxy formation scenarios. The halos consist of intermediate-age/old stars, 

I confirming the results of other studies that have shown the starburst phenomenon to be very centrally concentrated. 

There is no evidence for stellar-population age/metallicity gradients within the halos themselves. 

Key words, galaxies: evolution - galaxies: individual (NGC 1569) - galaxies: individual (NGC 4449) - galaxies: irregular 
- galaxies: photometry - galaxies: dwarf - galaxies: stellar content 
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^ ' 1. Introduction These local analogs of distant starburst galaxies provide the 

- - - opportunity for a very detailed study of galaxy formation 

The currently accepted theory of galaxy formation is that since we are able to resolve individual stars in these galax- 

galaxies form hierarchically, with the growth and evolution ies with the Hubble Space Telescope (HST) and thereby 

of large systems driven by build-u p through accret ion and build color-magnitude diagrams (CMDs), which are essen- 

merging of smaller units (see e.g. ICole et al.l[T99l . These tial tools for determining ages and abundances of the un- 

processes trigger intense, centrally concentrated bursts of derlying stellar populations. 

star formation (SF) that provide chemical enrichment and In the local universe, starbursts are typically found in 
thermal and mechanical heating of both the interstellar dwarf irregular galaxies. They are usually gas-rich systems 
and intergalactic medium (|Heckmanlll998t) . While starburst with a chemical comp osition comparable to that of distant 
galaxies are thought to be common at high redshift, the ob- primeval galaxies (e.g. lThuan fc Izotovlll999t ). Investigating 
servational information available for such distant galaxies the properties of underlying older populations in present- 
is very limited owing to their unresolved nature. In order day starbursting systems is the key to understanding their 
to gain better understanding of the processes that govern evolution over the cosmological timescale. The presence of 
galaxy formation we turn to nearby starbursting systems, older populations is a proof that a given galaxy has also 
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experienced episodes of SF in the past fe.g. lTosi et aLlflOQlL 
lGrebell[T998h . 

A number of recent studies have aimed at describing 
the properties and the origin of extended stellar halos that 
con tain these o l der populations in various galaxy types (see 
e.g. iSeth et al.l[2007L iSeth et al.ll200l Ide Jong et ahl [20091 
for the results on giant spiral galaxies from the GHOST 
project, Tikhonov | 2006a for the massive irregular M82, 
IVansevicius et al. 12004 for the very low-mass dirr Leo A). 
Recent advanced simulations allow one to ma ke predictions 
rega rding halo chemical abundance trends (jZolotov et al.l 
12010 ') and their surface brightness profiles following differ- 
ent formation scenarios (e.g. Cooper et al.l[201Q') : they also 
show tha t the halos c ould have formed both as a result of 
mergers (jBekkil I2008D or because of internal processes in 
isolated galaxies ( Stinson et al.l 120091 ). 

NGC 1569 and NGC 4449 are two of the closest and 
strongest starburst galaxies, at a distance of 3.04 Mpc 

I Grocholski et al. 2 011, in preparation) and 3.82 Mpc 
Annibali et al.l 120081 ) . respectively. Both have been stud- 
ied extensi vely w ith HST before (e.g. fc rcggio ct al. 1 9981 
lAloisi et al.ll2 00l'. A ngeretti et al.ll2005rfelair fc Fesenlll998l 
and lGelatt et al.ii20Qlr) and are known to display interest- 
ing and unusual properties. NGC 1569 has a star-formation 
rate (SFR) per unit area that is 2-3 times higher than 
in other strong starbursts and 2-3 orders o f magnitude 
higher than in Local Grou p dwarf irregulars (jAloisi et al.l 
l2001llMcQuinn et al.ll2010l ). Its three super star clusters are 
among the most massive and luminous in the local universe, 
making the galaxy an extreme example of a starburst. NGC 
4449 with the recent SFR of ~lMo/yr (jMcQuinn et al.l 
l201Clt) . is the only local "global" starburst, in that star for- 
mation aggears to occur throughout much of the galaxy 
(|Huntedll997l) . 

The recent anal ysis o f de ep HST/ACS photometry 
bv IGrocholski et al.l pOOj) and lAnnibah et al.l (,2008) (see 
Figures [Tj and [5| for mosaicked ACS images) revealed a va- 
riety of stellar ages and metallicities in both galaxies, in- 
cluding well-defined populations of red giant branch stars. 
With WFPC2 data that were obtained in parallel it is pos- 
sible to extend this analysis to much larger galactocentric 
radii, which is the topic of the present paper. We detect and 
quantify the presence of old stars at large radii (Section 2). 
We study stellar population gradients through variations in 
the Red Giant Branch (RGB) color and the relative num- 
ber of carbon stars as a function of radius (Section 3). An 
analysis of density profiles as traced by RGB star counts 
allows us to address the question whether the stars at large 
radii form a halo component or an extension of the inner 
disk (Section 4). Finally, we discuss the results and compare 
them with the literature (Section 5). 



2. Observations 

2.1. Fields 

We have observed both NGC 1569 and NGC 4449 with 
ACS and WFPC2. In both cases, ACS was centered on 
each galaxy with parallel WFPC2 imaging of a region ~ 6' 
away from the galactic centers. 

Our NGC 1569 data were taken during the HST Cycle 
15 program 10885 in November 2006 and February 2007 
(PI: A. Aloisi). Between the two sets of observations, HST 
was rotated roughly 90 degrees so that, while ACS was 



Fig. 1: ACS 3-color F658N (Ha), F814W (/) and F606W 
{V) composite image of NGC 1569 showing the locations 
of ACS subfields. Owing to the decrease in stellar density 
in the outer portions of the galaxy, we increased the size 
of our boxes as we move outward from the center of the 
galaxy. 
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Fig. 2: ACS 4-color F435W (B), F555W (V), F814W (/), 
and F658N (Ha) composite image of NGC 4449 showing 
the locations of ACS subfields. See the caption of Figure [1] 
for the explanation of different box sizes. 



centered on the galaxy, the parallel WFPC2 observations 
image two distinct regions in the outskirts of NGC 1569. 

NGC 4449 was imaged as a part of the Cycle 14 program 
10585 in November 2005 (PI: A. Aloisi). To fully image 
the galaxy's center, ACS data were taken at two adjacent, 
but overlapping, positions. As a result, the parallel WFPC2 
images are also adjacent to each other and overlap slightly. 

For both galaxies the WFPC2 images were taken in 
F606W and F814W broadband filters, which roughly cor- 
respond to standard V and / filters in the Johnson- Cousins 
system. The ACS images were taken in F606W/F814W 
for NGC 1569 and in F555W/F814W for NGC 4449. We 
here present the first photometric analysis of the WFPC2 
data and combine this with the previously published ACS 
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Fig. 3: Hubble Legacy Archive footprints: ACS and WFPC2 fields (black) for NGC 1569 (left) and NGC 4449 (right) 
pointings overlaid on Digitized Sky Survey (DSS) images. The fields depicted in light gray come from programs other 
than those listed in Table [T] and are not discussed in this paper because they do not go deep enough to contribute to our 
analysis. 



data for the ga lactic centers (jGrocholski et al.l 1200 jFI and 
lAnnibah et"al]l2008D . 

Table [T] lists the coordinates (a and S), adopted major 
axis position angles, axial ratios and distances, number of 
frames, total exposure times, galactocentric radii and el- 
liptical radii (the latter being the major axis radius of the 
corresponding isophote). The exact locations of all WFPC2 
fields are shown in Figure [3] in black. 

Other WFPC2 pointings are also available from the 
HST archive (programs 6111, 6253, 6423, 7909, 8059, 8133, 
8544 9244, 9249 and 9634 for NGC 1569 and 5446, 5971, 
6716, 8601, 9244 and 10522 for NGC 4449), as shown in 
Figure [3] However, only the data from the GO-8059 pro- 
gram (PI: S. Casertano) are deep enough for our purposes 
and are used here. 

2.2. Reduction 

For each filter and pointing the WFPC2 frames were pro- 
cessed through the HST data archive on-the-fiy repro- 
cessing system, using the most up-to-date reference files 
and software. They were then combined into a single mo- 
saicked image us i ng th e software package MultiDrizzle 
(jKoekemoer et all l2002|) to improve the sampling of the 
point-spread function (PSF) as well as to remove cosmic 
rays, bad pixels and correct for geometric distortions. In 
order to decide on the optimal pixfrac (the linear size of 
the "drop" in input pixels) and pixscale (the size of output 
pixels) values the so called "grid experiment" was carried 
out, following the recipe given in lFruchter fc SosevI (j2009tl . 
Its purpose was to determine what combinations of the two 
parameters provide the best resolution and PSF sampling 
for our dithered images. The general approach is that the 
pixfrac should be slightly larger than the pixscale value to 

^ For all calculations in the present paper we use an improved 
p hotometric catalog based on a reanalysis of the data presented 
in lGrocholski et al.l (|2008l ). 



allow parts of each "drop" to spill over to adjacent pixels. 
We found that for most of the images pixfrac— 0.7 and pixs- 
cale=0.6 provided the best resolution and PSF sampling. 
An example of a drizzled image is shown in Figure 21 




Fig. 4: Drizzled combination of 32 1000/llOOs images for 
NGC 1569, field 10885/1, F606W, WF3 chip. 



2.3. PSF photometry 

Point-source photometry was performed with the stand- 
alone version of the DAOPH OT/ALLSTAR stellar photom- 
et rv package (IStetsonlll987D following the method described 
in iGrocholski et al.l (j2008t) . A set of plots was created for 
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Notes. Adopted galaxy center coordinates (a and 5), major axis position angles (PAmaj), axial ratios (q), distances and corre- 
sponding arcmin/kpc scales; number of frames, total exposure times and galactocentric radii (r) and elliptical radii (m) for each 
WFPC2 pointing. The axial ratios and position angles were calculated by taking typical values from the IRAF ellipse task isophote 
fits at 2 arcmin radii. Elliptical radii m were calculated according to the following formula; = + [y^ /q^) where x and y are 
projected major and minor axis distances, respectively. 

References. (1) Grocholski et al. (2011, in preparation); (2) lAnnibali et al] ()2008f ). 



each pointing showing magnitude versus a (the uncertainty 
on the magnitude) , (the residual per degree of freedom of 
the PSF-fitting procedure) and sharpness (a measure of the 
intrinsic size of the object with respect to the PSF where 
indicates a "perfect" star, a negative value is indicative of, 
e.g. bad pixels or cosmic rays, and a positive value indicates 
extended objects such as faint background galaxies or unre- 
solved star clusters). An example is shown in Figure[5] where 
we can see that the a increases toward fainter magnitudes 
as expected. The values increase towards brighter magni- 
tudes, which is because systematic errors in the PSF model 
are more statistically significant at lower noise levels. The 
sharpness distribution widens toward fainter magnitudes, 
likely indicating increasing contamination from non-stellar 
sources. 

2.4. Aperture and CTE corrections 

The PSF magnitudes computed within ?alO pixel radius 
were corrected to the standard (X'S aperture for which the 
filter magnitude zeropoints are known. The second step 
was to correct these values to the hypothetical infinite 
aperture, which was done by applying a 0.1 mag offset 
(|Holtzman et al.lll995. ). 

Corrections for imperfect charge transfer efficiency 
CTE ) were calculated using the prescription by iDolphinl 
20091) . The computed corrections for most stars were mod- 
est, owing to the high sky backgrounds resulting from long 
exposures in broad-band filters. However, the corrections 
were as high as w0.14 mag for faint stars located far from 
the readout amplifier. 



1.2 7- 
1.0 - 




18 20 22 24 25 28 

mog 



Fig. 5: Distribution of the DAOPHOT parameters a (top), 
sharpness (middle), and (bottom) as a function of 
F814W VEGAMAG magnitude (a rough approximation of 
the Johnson-Cousins I) for NGC 1569, field 10885/1. 



2.5. Johnson-Cousins transformation 

The conversion of count rates DN/s into the standard 
Johnson-Cousins sys tem w as performed following the for- 
mulae by Holtzman et al.l ([1995) but using an updated 
set o f zeropoints and transformation coefficients (Dolphin 
l2009l) . The transformation formula is of the following form: 

SM = -2.5-log{DN/s)+Ti-SC+T2-SC'^ + Z+2.5-log{GR),{l) 



where SM is the standard magnitude, SC the standard 
color, Tl, T2 the transformation coefficients, Z the zero 
point and GR the gain ratio (slightly different for each 
chip). It was also necessary to specify multiple fits for both 
filters to accurately cover the entire color range: the trans- 
formation formula uses different coefficients for the color 
ranges <2.0 and >2.0. 
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2.6. Color-magnitude diagrams 



Figures [5] and [7] show CMDs of two of our ACS subfields 
(innermost and outermost) as well as all of the WFPC2 
fields for NGC 1569 and NGC 4449, respectively. For the 
WFPC2 fields we used the photometric quality parame- 
ters provided by DAOPHOT to remove false detections. 
Because the fields are sparse, fairly lenient cut values were 
chosen to maintain clear CMDs. 

In the case of NGC 1569 we applied the following cuts 
to each WFPC2 catalog: a < 0.3 and —1.0 < sfiarpness < 
1.0. No cuts in were applied since this would reject many 
real stars at the bright end of the CMD. A final check was 
made by visually inspecting the frames with "cleaned" cat- 
alogs overplotted to make sure that no spurious objects 
were present or no real stars were being rejected. The cata- 
logs contain 1176, 136, 675 and 391 stars for fields 10885/1, 
10885/2, 8059/1 and 8059/2, respectively. 

For the NGC 15 69 ACS subfields w-e use d the same pho- 
tometric cuts as in iGrocholski et al.l ()2008f ). The subfields 
contain 7570 (Al), 11173 (A2), 8570 (A3), 3441 (A4) and 
926 (A5) stars. 

For NGC 4449 there was little evidence for bright stars 
in the WFPC2 fields. This allowed us to adopt more strin- 
gent cut values than for NGC 1569. Applying the cuts 
a < 0.4, —0.5 < sfiarpness < 0.5 and < 3-0 allowed 
us to make sure that we were neither cutting out real stars 
nor keeping a lot of false detections. The final cleaned cata- 
logs contain 827 and 1681 stars in the first and second field, 
respectively. 

The four NGC 4449 ACS subfields contain 4043 (Al), 
11092 (A2), 10897 (A3) a nd 4510 (A4) stars. Again, the 
original cuts were applied (jAnnibali et al.|[2008t ). 
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2.7. Completeness 

The completeness for all WFPC2 and ACS regions was es- 
timated by performing artificial star experiments: artificial 
stars were added to the frames, each with a randomly as- 
signed magnitude from the range observed in the real data. 
Then, the same photometry steps were taken as for the 
real data. The completeness was calculated in the individ- 
ual filters, and then the final completeness was obtained by 
multiplying the completeness factors in the different filters 
(since in our photometry we require a star to be detected 
independently in each filter) . The number of recovered stars 
divided by the number of stars used as the input provided 
the estimate of completeness as a function of magnitude 
and position. Both the ACS and WFPC2 50% complete- 
ness limits are shown in Figures El and [71 



2.8. ACS photometry 

Here we provide a short summary of the ACS photometric 
reduction for NGC 1569 and NGC 4449. Spe c ific d etails 
can be foun d in Section 2 of Grocholski et al.l (|2008[ ) and 
Section 2 of lAnnibah et al.l (|2008[ ). 

Total exposure times for NGC 1569 were 61716 s in 
the F606W fiher and 24088 s in the F814W fiher. For 
NGC 4449, the total integration times in F435W, F555W, 
and F814W were ~ 3600 s, 2400 s, and 2000 s, respectively. 
The individual images were processed through the ACS 
pipeline (CALACS) and combined into final images us- 



Fig. 6: CMDs of the innermost (top left) and outermost 
ACS (top right) subfield and the four WFPC2 fields (mid- 
dle and bottom) for NGC 1569. The dashed-dotted lines 
indicate the calculated TRGB values and the long-dashed 
lines show 50% completeness levels. The two boxes shown 
in the 10885/1 plot indicate adopted color and /-magnitude 
ranges for selecting RGB and carbon stars (see Sections 13.31 
13.41 and 14.11 for details) . Smaller plot symbols are used for 
the Al subfield to avoid its saturation and make the differ- 
ent CMD features visible. 



ing th e MULTIDRIZZLE software package l|Fruchter et al.l 
120091) . 

Using the DAOPHOT/ ALLSTAR software (jStetsonl 

Il987h . we performed photometry in the following manner. 
The PSF template was computed from the most isolated 
stars in the frame to accurately model any spatial varia- 
tions in the shape of the PSF. We then used ALLSTAR 
to fit the PSF to independent source detection lists for the 
images in all bands. The photometry lists were then cross- 
correlated, and the resulting catalogs were cleaned of false 
detections, background galaxies, or stars with bad measure- 
ments by applying both positional and photometric qual- 
ity cuts. Stars that fell in the gap between the two ACS 
chips or near the edges of the array were removed because 
dithering caused these areas to have lower exposure time 
and therefore more noise. Detections near bright stars were 
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Fig. 7: CMDs of the innermost (top left) and outermost 
ACS (top right) subfield and the two WFPC2 fields (bot- 
tom) for NGC 4449. See the caption of Figure E] for further 
details. 



also removed because diffraction spikes tend to cause false 
detections. 

We corrected our photometry for C TE losses following 
the prescription of lRiess fc Mackl (|2004f ) and then converted 
instrumental magnitudes into the Jo hnson- Cousins system 
following the procedure outlined in ISirianni et al.l (|2005l) . 
Our final catalogs contain 369,213 objects with a measured 
magnitude in both V and I in the case of NGC 1569 and 
402,136 stars in the case of NGC 4449. 



3. Stellar population gradients 

3.1. CMD morphology and foreground contamination 

The innermost ACS CMDs (top left panel in Figures [5] and 
[7]) show all the evolutionary features expected at the mag- 
nitudes sampled by our data: a blue plume consisting of MS 
stars and blue loop stars at the blue end of their loops, a red 
plume consisting of red supergiants and oxygen-rich asymp- 
totic giant branch (AGB) stars, blue loop stars on tracks 
that cross between the two plumes, a horizontal finger of 
carbon stars extending redward from the tip of the RGB 
(TRGB), and a prominent RGB extending from its tip to 
the faintest detectable magnitudes. The first three features 
are due primarily to young stars (< 10^'^ yr), whereas the 
latter two are primarily due to intermediate- age and old 
stars (> 10^-5 yr). 

The starbursts in NGC 1569 and NGC 4449 are con- 
centrated toward the galaxy centers. The outermost ACS 
field of NGC 1569 (10885/A5) shows evidence only for the 
oldest CMD features of the RGB and carbon stars. These 
are also the dominant features in the outermost ACS field 




Fig. 8: NGC 1569 CMD for field 10885/1 (left) and the 
Besancon model predictions for foreground population for 
the same field (right). The boxes indicate a region popu- 
lated exclusively by foreground stars that was used to cal- 
culate the foreground scaling factor fsee 13. 4p . 



of NGC 4449 (10585/A4). However, that field stiU shows a 
small blue plume due to younger stars at V — I fa 0. This 
is the outer tail of the younger stars in the main body of 
NGC 4449, as quantified by the analysis of the spatial de- 
pendence of the stellar populations in Annibali et al. (2008; 
their figures 8 and 19). No evidence of young CMD fea- 
tures is present in any of the large-radii WFPC2 fields that 
form the primary focus of the present study. We have not 
attempted to set rigorous limits on any small amounts of 
young stars at large radii, but no such stars are needed to 
explain the features evident in the data. 

The NGC 1569 WFPC2 CMDs show a red giant branch 
of stars with colors {V-I) in the «1.5-2.5 range and a tip 
at / «24.5, a red tail of carbon stars extending above and 
to the right of the TRGB (2.4 < V-I < 4.0, 23.4 < / < 
24.4) mixed with a strip of foreground stars located between 
approximately V-I— 1.5, 1=18 and V-I^A, 1=25. In the 
case of NGC 4449, the red giant branch of stars has colors 
{V-I) in the 0.75-1.75 range and a tip at / w24.0 and the 
carbon stars are located approximately in the 1.75 < V- 
I < 3.0, 23.0 < / < 24.0 region. 

NGC 1569 is located close to the Galactic plane {b = 
11.24°), therefore our fields suffer from significant fore- 
ground contamination. The contrast of this foreground pop- 
ulation with respect to the stars in the target galaxy in- 
creases as one moves further out in the halo. To quan- 
titatively model the foreground population we used the 
Besancon models o f stellar populati on synthesis of the 
Milky Way Galaxy (jRobin et al.ll2003D . We generated syn- 
thetic CMDs for each of the fields to see how well they 
match our observational data (see example in Figure [5]). 
The simulation results fit the observational CMDs quite 
well when it comes to the location of foreground stars on 
the CMDs but they fail to agree in the exact number of 
counts. The predicted number of stars is ~30% higher than 
the number of foreground stars we observe in each field. 
This is, however, to be expected as the model is less ac- 
curate for low galactic latitudes and does not account for 
highly or irregularly absorbing clouds. 

For NGC 4449 (6 = 72.40°) the Galactic foreground 
contamination is negligible. The Besancon model predicts 
about six stars per square arcminute down to 7=28 mag. 
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3.2. TRGB magnitude 

The absolute /-band magnitude of the TRGB is ap- 
proximately a standard candle, with M j = —4.0 ±0.1 
(|Bellazzini et al.ll2004 [Barker et all2004D . For ages >2 Gyr 
and metallicities (Z) < 0.25 Zq there is little variation in 
the Mj with age or Z. Therefore, the observed TRGB mag- 
nitude is primarily an indicator of distance, with reddening 
also playing a role. 

Grocholski et al. (2011, in preparation) and 

lAnnibaU et"all (|2008D found the TRGB apparent I- 
band magnitudes for NGC 1569 and NGC 4449 to be 
24.47 ± 0.04 and 24.02 ± 0.05, respectively. Following the 
approach of the previous authors we calculate the TRGB 
magnitude for each of our WFPC2 fields using the method 
de veloped by R. P. van der Marel and described in detail 
in lCioni et all (|2000[ ). For N GC 4449, both o f our fields ex- 
cellently agreement with the lAnnibali et ahl ()2008il results, 
with an average of Itrgb = 24.06 ± 0.07. In NGC 1569, we 
find some variation between the TRGB brightness in each 
of our four fields, which could be caused by a difference in 
distance or reddening amongst the fields. At the distance 
of NGC 1569, our outermost field has a projected distance 
of 6.35 kpc. If we assume a thin disk geometry that is 
inclined 45° with respect to the line-of-sight, the difference 
in brightness between the TRGB in the outermost field 
and in the galactic center is ^0.005 mag. This is well below 
our photometric detection limits, therefore suggesting that 
reddening may be responsible for differences in the TRGB 
brightnesses. 

The foreground extinction maps of lSchleeel et al.l (|1998[ ) 
have a resolution of 6.1 arcmin and therefore provide insight 
into variations in foreground extinction between our fields. 
In general, we find little variation within each galaxy, with 
E(V-I) varying by less than ±0.02 mag for three of the NGC 
1569 fields and by less than ± 0.001 mag for both of the 
NGC 4449 fields. This is well within the error bars of the 
measurements and can therefore be safely ignored. The only 
exception is field 10885/1, which has an E(V-I) that is 0.15 
mag higher than the other NGC 1569 fields. We therefore 
correct for this differential reddening in all measurements 
of the 10885/1 field throughout the remainder of this paper 
to enable a fair comparison to the results of the other fields. 

When we correct our TRGB measurement in the 
10885/1 field for the differential reddening, we find the av- 
erage of aU four NGC 1569 fields is Itrgb = 24.47 ± 0.02, 
which excellently agrees with previous results. This agree- 
ment provides a useful consistency check on the photomet- 
ric calibration, RGB identification, and reddening estimates 
for our WFPC2 fields. The individual TRGB determina- 
tions are shown as horizontal dot-dashed lines in Figures [5] 
and [3 

3.3. RGB color 

In order to explore stellar population properties as a func- 
tion of distance from the galaxy center we used the WFPC2 
fields together with the ACS subfields shown in Figures [1] 
and [2 The WFPC2 fields for NGC 4449 lie on the same 
side of the galaxy, so we chose the ACS subfields along the 
diagonal of the ACS image that points roughly toward the 
WFPC2 fields. The WFPC2 fields for NGC 1569 are dis- 
tributed azimuthally around the galaxy, and we chose the 



ACS subfields somewhat arbitrarily to lie along the direc- 
tion toward the top left of the ACS image. 
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Fig. 9: RGB color variation as a function of elliptical radius 
TO from the galaxy center. ACS data points are shown as 
diamonds and WFPC2 data points as asterisks. The mea- 
sured colors shown here are not corrected for foreground or 
(possible) internal extinction. Best fits are shown as solid 
lines and their error ranges are indicated by long-dashed 
lines. 



To estimate the RGB color, boxes of stars with 
1.60 < V-I < 2.40, 24.4 < I < 25.4 (NGC 1569) and 
0.75 < V-I < 1.75, 24.0 < I < 25.0 (NGC 4449) were ex- 
tracted (see Figures E] and [T]). The stars were then binned 
according to color and a Gaussian was fitted to the color 
distribution. Its peak was taken as an estimate of the aver- 
age RGB color and its width divided by the square root of 
the number of stars as an estimate of the error. All RGB 
counts were corrected for the incompleteness determined as 
in Scction ETTl The results are shown as function of elliptical 
radius in Figure [SI 

Near the center of the galaxies the RGB colors are diffi- 
cult to interpret for two reasons. First, the presence of blue 
loop stars contaminates the blue side of the RGB, which 
leads to a color that might be biased on the blue side of the 
true color (as discussed in lAnnibali et al.l [20081 . Second, 
there is significant and variable extinction present intrinsic 
to the target galaxy, as evidenced by patch y dust obscura- 
tion (see e.g. iKobulnickv &: Skillmanlll997() . This may lead 
to an RGB color that is biased on the red side of the true 
color. Specifically, the central dip in the NGC 4449 profile 
is likely caused by the first effect. For these reasons we do 
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Table 2: Age and metallicity gradients. 





NGC 1569 


NGC 4449 


age 


3 Gyr 


10 Gyr 


3 Gyr 


10 Gyr 




0.00281±0.00031 


0.00137±0.00015 


0.00370±0.00030 


0.00178±0.00040 


Zgrad(Z/kpc) 


0.00007±0.00007 


0.00002±0.00004 


0.00005±0.00033 


0.00002±0.00015 


agegrad(Gyr/kpc) 


0.16±0.19 


0.45±0.64 


0.11±0.45 


0.50±1.50 



Notes. Metallicity, metallicity gradient (assuming ages of 3 Gyr and 10 Gyr, respectively) and age gradient (assuming constant 
Z) estimates for each galaxy, based on the observed RGB color (gradient). 



not interpret the RGB colors in the central 2.5 arcmin of 
each galaxy. 

For both galaxies, the color stays relatively constant but 
with some up and down variations that exceed the uncer- 
tainties in the photometry. To quantify this we performed 
a straight-line fit to the data outside of 2.5 arcmin. We cast 
the fit in the form c = A + B ■ (to— < m >), where < to > is 
the average m of the fields (both in kpc), c is the color, and 
A and B are the average RGB color and color gradient per 
kpc, respectively. There are usually two methods to set the 
weights in linear fits, one based on the known random un- 
certainties in the data points, and the other based on the 
rms of the points around the best fit. In our analysis we 
used the higher of these two values. This ensures that any 
unknown systematic uncertainties are accounted for in the 
error bars we derive on A and B. The best fits together with 
their error ranges are shown in Figure For both galax- 
ies the data are consistent with zero RGB color gradient 
outside 2.5 arcmin, with the maximum allowed gradients 
being -0.002 or -h0.012 mag/kpc for NGC 1569, and -0.017 
or -1-0.025 mag/kpc for NGC 4449. 



The RGB color is a highly degenerate indicator of stel- 
lar population properties, depending on extinction, metal- 
licity, and age. The goal of the present analysis is, there- 
fore, to search for relative population gradients and not 
absolute population properties. Nonetheless, it is useful to 
briefiy asses the populations that may be consistent with 
the observed color. For this we use the stella r popu lation 
synthesis model predictions of iGirardi et al.l (|2002[ ). For 
NGC 156 9 we use the for eground extinction E{V — /) = 
0.60 from Burstein fc Heiles (.1984 ) (see the discussion in 
iGrocholski et al.ll2008ll and for NGC 4449 w e use the fore- 
roun d extinction E{V — I) = 0.019 from ISchlegel et al.l 
19981) . We assume that the internal galaxy extinction is 
negligible beyond 2.5 arcmin. 



Table [5] shows what the measured RGB color implies for 
the metallicity Z if one assumes either age = 3 Gyr or 10 
Gyr. This can be compared to the metallicities of the gas in 
the central region, which are Z »0.005 and 0.0 06 for NGC 
1569 and NGC 4449 (see, e.g. iGreggio et all (^1998) and 
lAnnibali et al.l (|2008[ ) and references therein). Regardless 
of the assumed age, the stellar metallicities found from an- 
alyzing the RGB colors are lower than the gas metallicities 
in each galaxy. This was to be expected, because the gas has 
likely been metal-enriched through recent star formation. 

Using these results we estimated Z gradients at fixed 
age and age gradients at fixed Z (see again Table [2]). The 
large error bars mean that the results provide only an upper 
limit to any gradient that may be present intrinsically. 



3.4. Carbon star counts 

Another CMD feature that allows the study of radial gra- 
dients in old populations is the red tail of carbon stars. A 
unique identification of carbon stars requires spectra but as 
an approximation we can select stars from the CMD using 
the region that is belie ved to be predomin antly composed 
of carbon stars (see e.g. lMarigo et al.ll2003[) . 
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Fig. 10: C/RGB stars ratio as a function of elliptical radius 
m from the galaxy center. ACS data points are shown as di- 
amonds and WFPC2 data points as asterisks. Best straight 
line fits are shown as solid lines and their error ranges are 
indicated by long-dashed lines. 



To obtain the estimate of the number of carbon stars 
in NGC 1569 we first determined a foreground scaling fac- 
tor by counting stars in a CMD region that contains only 
foreground stars (see Figure [S]) in both the data and the 
Besancon model and dividing the former by the latter. 
Then, we subtracted from the data the foreground model 
predictions multiplied by this scaling factor. 
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Carbon stars were selected at 2.40 < V-I < 3.85, 
23.4 < / < 24.4 (NGC 1569) and 1.75 < V-I < 3.0, 
23.0 < / < 24.0 (NGC 4449). The color limit was chosen 
to avoid a s ignificant contributio n of RGB and oxygen-rich 
AGB stars (jAnnibali et al.|[2008f) . The chosen ranges are de- 
picted as boxes in the middle left panel of Figure [S] and the 
bottom left panel of Figure [T] Figure [TU] shows the results 
as a function of elliptical radius. The same completeness 
corrections were applied as for the RGB color. For simplic- 
ity, we refer to the ratio of the counts in the CMD boxes 
as the C/RGB ratio. However, it should be kept in mind 
that no attempt was made to to uniquely classify the stars 
in the chosen boxes. 

Using the C/RGB ratio as an absolute indicator of stel- 
lar population properties is difficult for several reasons. For 
example, the number of carbon stars depends on the life- 
time of the thermally-pulsating (TP) AGB phase, which is 
only poorly known. Also, many carbon stars c ould be dust 
enshro uded and unseen in the /-band; see e.g. I B over et al.l 
(|2009t) , who have revealed that up to 40% of AGB stars may 
be missing from optical CMDs. However, the C/RGB ratio 
can again be used as an indicator of relative population gra- 
dients. In any stellar population model the number of TP- 
AGB and carbon stars tends to decrease compared to the 
number of RGB stars with increasing population age (e.g. 
[Thuan fc Izotovl [2005). Therefore, gradients in the C/RGB 
ratio might be indicative of age gradients. 

As before, we do not interpret the results inside 2.5 ar- 
cmin, because the "RGB box" in the CMD may be contam- 
inated by younger stars that are not actually on the RGB. 
Outside of the central 2.5 arcmin the observed gradients 
are very shallow. The straight line fits in the figures show 
that both galaxies are consistent with zero gradient. 



4. Density profiles 

4.1. RGB counts 

Observational evidence indicates that different types of 
systems follow different density profiles. For example, the 
metal-weak halo of the Milky Way follows a roug hly power- 
law density profile: p ~ r'^-S (jJuric et al.l 120081) . Elliptical 
galaxies tend to have profiles to large radii that resemble 
the d e Vaucouleurs law : surface brightness is linear with 
r^^'^ (jSchade et al.l[T999^ . Spiral galaxy disks tend to have 
an exponential fall-off with radius: their surface brightness 
is linear wi th r (and they often have a truncation at some 
radius; e.g. lde Jong et al.ll2007l ). These different profiles are 
straight lines in plots of log density versus either log{m), 
m^^^ or TO, where m is the elliptical radius. 

Because we are interested in understanding to what kind 
of population our stars belong, we study their density pro- 
files using the number density of RGB stars at different 
radii. To do this, we use the same boxes and completeness 
corrections as we did to study the RGB color. In Figure [TT] 
we plot the density profiles of NGC 1569 (left) and NGC 
4449 (right) as a function of the elliptical radius, to (top), 
as well as m^^^ (middle) and log(m) (bottom). The equa- 
tion for the best straight line fits are given, with the fit 
overplotted on the data. For the exponential fits, the disk 
scale lengths, Rd, are given and Re, the effective radii, are 
indicated for the to^/'' fits. 

The values of Rd quoted in the figures are higher 
than the ones previously given in the literature (0.385' for 



NGC 1569 (|\yillett et al.ll2005l) and 0.763' for NGC 4449 
([Hunter et al.l loO^Tboth converted from values in kpc us- 
ing the distances adopted by the authors) . This is presum- 
ably because previous authors measured integrated light in 
the optical, where one is sensitive to young stars that are 
more concentrated toward the galaxy center than the RGB 
stars studied here. 

For both galaxies the star counts appear to be fitted rea- 
sonably well with the de Vaucouleurs and power-law profiles 
(middle and bottom panel of Figure [TTl respectively). The 
log(7V) vs. TO plots (Figure [TTl top panel) show that the 
number of stars at large galactocentric distances is higher 
than a simple exponential model would predict. This sug- 
gests that the stars we see at large radii are not just the 
extrapolation of an inner exponential disk. Instead, they 
are probably part of a spheroidal stellar halo component. 

While in principle we are able to fit a double exponential 
to the data (see Figure fT2)) . this does not mean that a thick 
disk is the favored explanation of the elevated density we 
see at large radii. Many of our fields are well away from 
the projected major axis, and in view of this a rounder 
component, i.e. a halo, seems more plausible. 

4.2. Integrated light profiles 

In order to perform a consistency check, we determined the 
profiles of integrated light from 2MASS i^T-band images for 
comparison to our RGB counts. K is the reddest band in 
the 2MASS survey, and is most dominated by light from 
RGB stars. Because intensity is linearly proportional to the 
number of stars, one would expect the RGB count profiles in 
Figure[Tl]to match the measured 2MASS surface brightness 
profiles after applying a vertical shift. 

The intensities at different major axis radii were cal- 
culated from the 2MASS images using the /i?^/[3 task el- 
lipse, with the ellipticity, PA and axial ratio values fixed 
to values typical at 2 arcmin radii (see Table 1). The cal- 
ibrated 2MASS magnitudes were then obtained by apply- 
ing the following formula: mag = MAGZP — 2.5 log]^o('5'); 
where the S is the integrated, background-subtracted flux 
in ''DN^^ (data numbers) and MAGZP is the zeropoint 
obtained from the image header. We truncated the profiles 
at ^0.4 % of the sky value, the same as the value adopted 
for profiles published in the 2MASS Large Galaxy Atlas at 
http://irsa. ipac. caltech. edu/. 

As can be seen in Figure [131 the resulting profiles line 
up continuously with the shifted RGB-count profiles, which 
was to be expected since i^-band measures mostly light 
from RGB stars. This also confirms that individual star 
counts are an excellent way of following density profiles to 
large radii, where the integrated light falls well below the 
sky background and becomes hard to quantify. 



5. Discussion and conclusions 

We obtained and analyzed HST/WFPC2 images of several 
fields in the outer halos of the nearby starburst galaxies 
NGC 1569 and NGC 4449. The fields were imaged using 
the F606W and F814W filters. We reduced the data and 



^ IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of 
Universities for Research in Astronomy (AURA) under coop- 
erative agreement with the National Science Foundation. 
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log(m) (orcmin) log(m) (orcmin) 

Fig. 11: NGC 1569 (left) and NGC 4449 (right) density profiles as a function of m (top), m^^^ (middle) and log(m) 
(bottom), where m is the elliptical radius from the galaxy center. N is the number of stars in the RGB boxes shown in 
Figures [S] and [T] The fits fixed to the ACS data points (top) and all points (middle and bottom) are shown as solid lines. 
The symbol types are the same as in Figures ^ and [TUl The dotted lines show disk scale lengths Rd (upper panels) and 
effective radii Re (middle panels). 



performed point source photometry to obtain CMDs of / 
vs. V — I. We used these CMDs to perform a study of stellar 
population and density gradients in these galaxies. To max- 
imize the radial range available for analysis, we included in 
our study also the CMDs (in the same bands) derived from 
our prior HST/ACS photometry near the galaxy centers. 
Corrections for foreground contamination, incompleteness, 
and differential dust extinction were calculated and applied 
as necessary. 

The HST/WFPC2 fields in our study are located 
roughly between 5 and 8 scale radii (exponential disk scale 
radius or effective radius) from the galaxy centers. The cor- 
responding elliptical (isophotal or major axis) radii fall be- 
tween 7 and 14 scale radii. Despite these significant radii, 
we detect stars belonging to these galaxies in all the fields. 
This implies that the galaxies have faint outer stellar en- 
velopes, and are not tidally truncated within the range of 
radii addressed by our study. Similar results have been 



obtained for other st ar-formin g dwa r fs, see, e.g. the re- 



3.g . the re- 

cent work by iNoel fc Gallart. ()2007D . iTikhonovl (|2006b[) . 
lAlonso-Garcia et all ^200^ . 



The CMDs of the HST/WFPC2 fields imply that the 
stars fall predominantly in the evolutionary sequences of 
the RGB and the TP-AGB (carbon star) phases. These 
phases correspond to stars of intermediate and old ages (in 
excess of a few hundred Myr). Although younger stars on 
the main sequence, blue loop, and red supergiant phases 
are present in abundance near the centers of these star- 
burst galaxies, they are not seen in the outer HST/WFPC2 
fields. This confirms results from our o wn previous work 
on the central regions of these galaxies (jGrocholski et al] 
l2008l lAnnibali et aLll2008D , as well as from studies of other 
galaxies in the literature (e.g. lTosi et al.l[2"00lLlTolstov et al.l 
2009) , that the starburst phenomenon tends to be strongly 
concentrated toward the galaxy center. The presence of 
older stars at all radii indicates that the starbursts occurred 
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in pre-existing galaxies with prior star formation, and were 
not the first star-formation episodes in an otherwise pri- 
mordial gas cloud. 

We identified and quantified the TRGB magnitude 
in all HST/WFPC2 fields. The average WFPC2 TRGB 
magnitudes (not corrected for foreground extinction) are 
24.47 ± 0.02 (NGC 1569) and 24.06 ± 0.07 (NGC 4449). 
These values excellently agree with the results obtained 
from our pr eviously analyzed HST / ACS data: 24.47 ± 0.04 
(NGC 1569;'Grocholski et aL|[200l and 24.02 ±0.05 (NGC 
4449; [Annibali et al. 20oCrespectively. The new TRGB 
results therefore support the galaxy distances derived in 
these previous studies. 

To study the origin and characteristics of the 
intermediate- age/old outer stellar envelopes of the sample 
galaxies (outside the central starburst), we measured the 
observed radial profiles of the RGB color and the number 
ratio of carbon to RGB stars. Gradients in these quantities 
would be tell-tale signs of radial gradients in either the age 
or metallicity of the population. However, to within the un- 
certainties of our analysis, no gradients are detected (the 
maximum allowed gradients are -0.002 or -1-0.012 mag/kpc 
for NGC 1569, and -0.017 or -^0.025 mag/kpc for NGC 
4449). This does not mean that no gradients exist, but 
merely that better data would be needed to detect any gra- 
dients that might exist. 

We measured the density profiles of RGB stars over the 
radial range accessible to our study. At small radii we show 
that the profiles are consistent with 2MASS measurements 
of near-IR integrated light. The profiles were fitted with 
exponential, de Vaucouleurs, and power-law models. While 
the latter profiles provide acceptable fits, exponential pro- 
files do not. Both galaxies have an excess of stars at large 
radii above what a simple exponential model would predict. 
This indicates that the intermediate- age/old stars at large 
radii at not merely the outer part of an exponential disk. 
Instead, the observed density profiles are more typical for 
an extended stellar halo. Such halos are believed to be a 
common byproduc t of hierarchical formation (but see also 
[St mson et al.ll2009l for a discussion on how such halos can be 
formed in isolated galaxies, i.e. without external perturba- 
tions). Stellar halos have been fou nd in all dwarf gala xies 
studied to a sufficient depth (see IStinson et al.l l2009l and 
the extensive reference list therein). Our study shows that 
NGC 1569 and NGC 4449 have similar halos to those seen 
in other galaxies of similar morphology as well as that deep 
resolved stellar photometry enables one to distinguish be- 
tween the profiles of the halos and the profiles of the star- 
forming disks they contain. 
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Fig. 12: Double exponential fits (solid lines) to the data, with their single exponential components (dashed lines) indicated 
(NGC 1569 on the left and NGC 4449 on the right). 
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Fig. 13: Comparison of surface brightness profiles from 2MASS K-band images with our RGB star counts for NGC 1569 
(left) and NGC 4449 (right). A vertical shift was applied to the RGB profiles to line the profiles up. 
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